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We report our findings on quantum phase transitions in cold bosonic atoms in a one dimensional 
optical lattice using the finite size density matrix renormalization group method in the framework 
of the extended Bose-Hubbard model. We consider wide ranges of values for the filling factors and 
the nearest neighbor interactions. At commensurate fillings, we obtain two different types of charge 
density wave phases and a Mott insulator phase. However, departure from commensurate fillings 
yield the exotic supersolid phase where both the crystalline and the superfiuid orders coexist. In 
addition, we obtain signatures for solitary waves and also superfluidity. 
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I. INTRODUCTION 

The supersolid phase, first reported in 4 He [l[ , is char- 
acterized by the coexistence of the superfiuid and crys- 
talline orders. This phase has beenpredicted in sev- 
eral bosonic lattice systems @, H, El IE IE 0, how- 
ever, there has been no unambiguous observation of this 
phase so far. Kim et al had reported its observation in 
solid 4 He M but a number of studies disagree with this 
claim [I 

In recent years, the advances in the manipulation of 
cold bosonic atoms in the optical lattices have opened 
up a new route to investigate quantum phase transi- 
tions [HI, EH • This approach has many advantages over 
the conventional solid state techniques, such as for exam- 
ple, flexibility in controlling the parameters and the di- 
mension of the lattice by tuning the laser intensity. A sys- 
tem of cold bosonic atoms in an optical lattice can be ad- 
equately described by the Bose-Hubbard model 0, HH . 
However, if the atoms possess long range interactions due 
to the presence of magnetic dipole moments, for exam- 
ple, then they could exhibit a number of different novel 
phases. In particular, the existence of such interactions 
could result in the supersolid phase @, 0, HI, E3| ■ The 
fairly recent observation of the Bose-Einstein condensa- 
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tion of 52 Cr atoms [18J, which have large magnetic dipole 
moments could ultimately lead to the observation of this 
unusual phase. 

In this context, we re-investigate the system of bosonic 
atoms with the long range interaction using the extended 
Bose-Hubbard model given by 

H = -t J2 (ojoj + H.c) 

<i,j> 

+ n ^ ni ~ 1 ) + v Uin i- w 

i <i,j> 

Here t is the hopping amplitude between nearest neigh- 
boring sites < i,j >. a\(a,i) is the bosonic creation (an- 
nihilation) operator obeying the Bosonic commutation 
relation [a], cij} = 5i.j and rii = ajaj is the number oper- 
ator. U and V are the onsite and the nearest neighbor 
interactions, respectively. We rescale in unit of the hop- 
ping amplitude, t, setting t = 1, making the Hamiltonian 
and other quantities dimcnsionless. 

In the absence of long range interactions, the model 
given in Eq. [I] reduces to the Bose-Hubbard model which 
exhibits a superfiuid (SF) to Mott insulator (MI) tran- 
sition at integer densities of bosons [l4j. However, for 
non-integer densities, the system remains in the super- 
fluid phase which is compressible and gapless. The Mott 
insulator phase, however, has finite gap and is incom- 
pressible. The extended Bose-Hubbard model, given in 
Eq. [TJ has been studied earlier using different meth- 
ods [I [H, [H [S| including DMRG ||3, [H H. The 
inclusion of the nearest neighbor interaction gives rise to 
the charge density wave (CDW) phase for integer and 
half integer densities @, [H, [H H H [H, H that has 
finite gap, finite CDW order parameter, vanishing com- 
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pressibility and a peak in the density structure function 
at momentum q = it. In the CDW phase the bosons 
occupy alternate sites and the unoccupied ones being 
empty. For example, when the density p = 1/2, the 
distribution of bosons has a | 101010 ■■•) structure 
while for p = 1 it is | 2 2 2 • • •). To distinguish 
between these two CDW ground states, the former is re- 
ferred to as CDW-I and the latter, CDW-II 0. This 
model has been studied recently using quantum montc- 
carlo resulting in the prediction of the supersolid phase 
when the density of the bosons is no longer commensu- 
rate. We re-investigate this model by departing from 
both half and integer filling for large and intermediate 
onsite interaction strengths and obtain the phase dia- 
gram using the finite size density matrix renormalization 
group(FS-DMRG) method [25j,l26| and throw more light 
on the supersolid and solitonic phases. 

This paper is organized as follows. In Sec. HH we will 
discuss the method of our calculation using FS-DMRG. 
The results with discussion are presented in Sec. [Hi] fol- 
lowed by our conclusions in Sec. IIVI 



gives information about the density distribution of dif- 
ferent phases and finally the existence of the CDW phase 
is confirmed by calculating its order parameter: 



II. METHOD OF CALCULATION 

To obtain the ground state wave function and the en- 
ergy for the system of N bosons on a lattice of length L, 
interacting via an on-site and a nearest neighbor interac- 
tion, we use the FS-DMRG method with open boundary 
conditions [25], [26| . This method is best suited for one 
dimensional problems and has been wi dely used to study 
the Bose-Hubbard model [H, H3, ^ . We have 
considered four bosonic states per site and the weights of 
the states neglected in the density matrix formed from 
the left or right blocks are less than 10 -6 (23J. In order to 
improve the convergence, at the end of each DMRG step, 
we use the finite-size sweeping procedure given in [23l. l25j| . 
Using the ground state wave function \iPln) and energy 
El(N), where N refers to the number of bosons and 
L, the length of the lattice, we calculate the following 
physical quantities and use them to identify the different 
phases. The chemical potential p of the system having 
density p = N / L is given by 



M = 



S E L (N) 
SN 



(2) 



and the gapped and gapless phases are distinguished from 
the behavior of p as a function of p [28| . The compress- 
ibility k, which is finite for the SF phase, is calculated 
using the relation 



Sp 

The on-site local number density (n^), defined as, 
(rn) = (ipLNlnilipLN), 



(3) 



(4) 



(5) 



When the ground state is a CDW, FS-DMRG calcula- 
tion with open boundary leads to an artificial node in the 
density distribution at the center due to left-right sym- 
metry. We circumvent this problem by working with odd 
number of sites. In our calculations, we start with five 
sites instead of usual choice of four sites and increase the 
length up to L = 101 adding two sites in each DMRG 
iteration (23|. After reaching the desired length L = 101, 
we vary the number of atoms N from 26 to 125 to scan 
a wide range of densities [28| . In this work we consider, 
two different values of the onsite interaction strengths: 
U = 5 and 10 and vary the nearest neighbor interaction 
strength V from to U. The choice of U is guided by an 
earlier work [23| where a direct MI to CDW-II transition 
for U = 10 and a MI to SF to CDW-II for U = 5 were 
observed as V is varied at a density p = 1. In this work, 
we extend this calculation to a wider range of densities 
and obtain a richer phase diagram consisting of super- 
solid and solitonic phases in addition to SF, CDW-I and 
CDW-II. We begin our discussion for U — 10 and later 
comment on our results for U = 5. 



III. RESULTS AND DISCUSSION 

It is well known that the Bose-Hubbard model (Eq. [1] 
with V = 0) has a superfluid ground state when den- 
sity p is not an integer and exhibits a quantum phase 
transition from the superfluid to the Mott insulator for 
integer densities [l4[ at a critical value of onsite interac- 
tion Uc that depends on p. (For example, Uc ~ 3.4 for 
p = 1 [13, HHO The Mott insulator has finite gap and 
zero compressibility while the superfluid is gapless and 
compressible. In the presence of a finite nearest neigh- 
bor interaction V, an additional insulator phase. CDW, 
appears at commensurate densities. As noted in |22l.[23|. 
a CDW-I occurs at p = 1/2 and at p = 1, depending on 
the value of V either a MI or a CDW-II appears. Since 
we are dealing with only an on-site and a nearest neigh- 
bor interaction, the commensurate densities for model in 
Eq. [T] are integers and half integers. We begin by study- 
ing the possible phases at commensurate densities, before 
we investigate the phases at incommensurate densities. 

The gapped phases are easily obtained from the depen- 
dence of the density p and the compressibility n on the 
chemical potential p. Figure [T] shows the dependence of p 
on p for a fixed value of U =10 and V ranging between 2 
and 10. The gapped phases appear as plateaus with the 
gap equal to the width of the plateau, i.e., p + —f/T , where 
p + and pT , respectively, are the values of the chemical 
potential at the upper and lower knee of the plateau. For 
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small values of V, Fig. Q] has only one plateau at p = 1. 
However, as we increases V an additional plateau ap- 
pears at p = 1/2. We calculate the compressibility using 
Eq. [3] and is also given as a function of p in Fig. [2] for 
two generic values of V. The smaller value, V = 2, has 
just one plateau at p = 1 while V" = 7, has two, at den- 
sities p = 1/2 and 1. As expected, the compressibility is 
zero over the range of p values where the plateaus occur, 
while it is finite elsewhere. The incompressible insulator 
and compressible superfluid regions can be separated out 
by picking up p + and pT and plotting them in the p, — V 
plane. From Figs. [T] and [H we see that, (i) a gapped 
phase occurs at p — 1/2 for V > 2.8, (ii) for p = 1, the 
gap remains finite for all values of V and (iii) the gap is 
zero for other values of p. 




FIG. 1: (Color online)The density p as a function of the chem- 
ical potential p for different values of V. The plateaus at the 
commensurate fillings indicate the existence of finite gap in 
the system. 

The nature of the compressible and incompressible 
phases can further be understood from the local den- 
sity distribution (n,) and the charge density wave order 
parameter Ocdw given by Eqs. [4] and [5] The variation 
of local density as a function of the lattice sites are given 
in Figs. [3] and O for densities around p — 1/2 and 1. At 
commensurate densities, say, p = 1/2, the charge density 
wave nature of the phase is clearly observed for V = 5.6 
in Fig. [2ICc) . Alternative variation of the density of bosons 
between one and zero is the signature of CDW-I phase, 
| 1 1 1 • • •) type 0,113. Similarly, for density p = 1, 
the density oscillation of the type | 2 2 2 • • •), for 
V = 9 suggest the CDW-II phase. From the gap, the 
compressibility and the density oscillations, we can con- 
clude that for U — 10 and p = 1/2, we have a SF to 
CDW-I phase transition at V ~ 2.8. However, for p = 1, 
there is no superfluid phase and the transition is from 
MI to CDW-II at the critical value V c ~ 5.4. These 




FIG. 2: (Color online) Variation of n and p with p for V = 2 
(top panel) and V = 7 (bottom panel). The plateau regions 
have zero compressibility while it is finite elsewhere. 



results are consistent with the earlier results in the liter- 
ature [USUI- 

We now turn our attention to the case when p is not 
commensurate to the lattice length, i.e., p ^ 1/2 or 1. 
From Fig. [51 we observe that the compressibility is al- 
ways finite for incommensurate densities indicating that 
these regions of the phase diagram correspond to the su- 
perfluid phase. However, the local density distribution 
and Ocdw reveal the richness of the phases present in 
the compressible regions of the phase diagram. Inter- 
esting phases appear when the nearest neighbor inter- 
action is large enough to obtain a CDW-I or CDW-II 
phase at commensurate densities. Let us first consider 
densities close to 1/2. When V is less than the critical 
value, Vc ~ 2.8, for the SF-CDW transition, we expect 
only the superfluid phase. However, for V > Vc the 
ground state shows solitonic behavior for densities close 
to p = 1/2. Fig. [3J shows the the local densities (rij) as a 
function of the lattice sites i for different densities. The 
panel labeled (c) corresponds to the commensurate den- 
sity p = 1/2 where we clearly observe the CDW nature 
of the ground state, (b) and (d) shows the density varia- 
tions of the ground state where one boson has been added 
and removed from the system at p — 1/2 respectively. 
Similarly, panels (a) and (e) show the density variations 
when two bosons have been added and removed respec- 
tively. The density profiles can be understood as follows. 
Moving away from commensurate densities, the solitons 
distort the periodic ground state by breaking the long 
range crystalline order as a modulation in the density 
wave that minimizes the ground state energy of the sys- 
tem dm. 

To understand the solitons we calculate the CDW or- 
der parameter for each unit cell. In contrast to the su- 



pcrfluid and Mott insulator phases that have one site per 
unit, the unit cell of CDW phase consists of two lattice 
sites. Referring to these two sites as 1 and 2, we define 
the CDW order parameter per unit cell as 



(panel (b)) at a density p = 1. 



O c cdw = (»i> - <na>. 



(6) 



The CDW phase has two degenerate ground states 
corresponding to the two local density distributions | 
1 1 1 • ■ ■ } and | 1 1 1 • ■ • ). The CDW 
order parameter, Ocdw> f° r these two degenerate states 
is equal to 1 and —1 respectively. Figure H] shows the 
Oqdw f° r tne same set 01 densities considered in Fig. [3] 
In Fig. [U the center panel (c) has density p — 1/2, while 
(b) and (d) represent the system in panel (c) with one 
boson added and removed respectively and (a) and (e) 
have two bosons added and removed with respect to (c). 
We notice that the Ocdw 1S uniform and close to one for 
p = 1/2. Since we work with odd number of sites with 
open boundaries, energy consideration leads to a CDW 
ground state which represents | 10 10 10 • • • } state. 
When we add or remove one boson from this state, we get 
two solitons that modulate the density distribution and 
break the long-range crystalline order. The extra particle 
or hole splits into two solitons of equal mass [29|, [3(j ■ The 
two solitons can move across the lattice without causing 
any energy, however, if we want to get rid of them, we 
need to spend lots of energy to flip the bosons. 

Similarly, removal or addition of two bosons result in 
four solitons. Continuing this process results in more soli- 
tons, until a critical density is reached, when the density 
oscillation completely dies out and the superfluid phase 
is obtained. Therefore, starting with the CDW-I phase 
and changing the density from its commensurate value 
of p = 1/2 by either adding or removing bosons, leads 
to solitons+SF phase that finally becomes a superfluid. 
The transition from solitonic to superfluid is more like a 
crossover rather than a real phase transition. The soli- 
tonic phase is obtained only very close to p = 1/2 and 
remains stable for the entire range of V considered on 
the hole side (p < 1/2). However, on the particle side 
(p > 1/2), the solitonic phase remains stable only up 
to some critical value of V — Vc ~ 6.4. For V > 6.4, 
doping below half-filling breaks the CDW ground state 
into a solitonic state that eventually goes into a super- 
fluid phase. However, this does not happen when bosons 
are added above half-filling. For example, the variation of 
(rii) as a function of sites i for three different densities are 
given in Fig. [5l for V = 9. The panel (b) represents the 
CDW-I phase at p — 1/2 while (a) and (c) respectively, 
correspond to the ground state with density obtained by 
removing and adding one boson to the CDW-I ground 
state. While a solitonic phase appears when bosons are 
removed (p < 1/2), Fig. [S^c) suggest that the CDW-I 
phase is robust for p > 1/2. Similar behavior is also seen 
when doping around p = 1. Fig. [5] shows (rii) as a func- 
tion of i for densities around p = 1. Panels (a) and (c) 
correspond to the density of the ground state obtained 
by removing and adding one boson to the CDW-II state 



0.5 



0.5 
0.5 
0.5 



(a) 



(b) 



(c) 



(d) 




FIG. 3: (Color online)Local density (n,) as a function of lat- 
tice sites i. Panels (a) and (b) shows the solitonic signature 
when p < 1/2. (c) shows the signature of CDW-I where every 
alternate site is occupied with one boson for p = 1/2. Panels 
(d) and (e) shows modulation of the CDW-I phase for p > 1/2 
and are once again in the solitonic phase. 
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FIG. 4: (Color online)CDW order paramater per unit cell, 
Ocdw as a function of sites, (a) and (b) exhibits the solitonic 
signature for densities p < 1/2. In panel (c), O c e D ' w is flat, 
thereby indicating CDW-I with each alternate site occupied 
with one boson, (d) and (e) show the signature for the solitons 
at densities p > 1/2. 

It turns out that for large V, the region between 
1/2 < p < 1, i.e., between CDW-I and CDW-II, always 
remains in the CDW phase even though the density is 
not commensurate to the lattice length. The CDW order 
in the system is determined by calculating the CDW or- 
der parameter given by Eq. [5] and is given in Fig. [7] For 
small values of V, the Ocdw is zero for all the densities 
except at p — 1/2 signalling the CDW-I phase. However, 
as V increases, an additional peak develops at p = 1 for 
V > 5.4 which corresponds to the CDW-II phase. The 



most interesting feature is that the Ocdw remains fi- 
nite in the region 1/2 < p < 1 for large values of V. It 
may be noted that the compressibility for 1/2 < p < 1 
is always finite. So the bosons move freely on the CDW 
background and prefer to occupy sites which are already 
occupied. Consider a system with one extra boson added 
to the ground state corresponds to CDW-I. If the added 
boson occupies an empty site, the energy cost is only due 
to the nearest neighbor interaction and is of the order of 
2V. However, if the added boson occupies a site which 
is already occupied by an another boson, the energy cost 
is due to the onsite interaction and is of the order of 
U which is relatively smaller than 2V for large V. The 
extra bosons therefore move between the occupied sites 
with a finite hopping amplitude leading to a long range 
correlation in the lattice, which yields to finite compress- 
ibility in the region l/2<p<las shown in Fig. [5] 
Similar behavior persists for doping above p > 1 as given 
in Fig. [5J 





FIG. 5: (Color online) Variation of the local density (nj) as 
a function of i for V = 9. Panel (a) shows the ground state 
density when a boson is removed from the state at p — 1/2 
shown in panel (b), which is a CDW-I and panel (c) corre- 
sponds to the state obtained by adding a boson to p = 1/2 
state. 

Therefore we can conclude that for small V apart from 
commensurate fillings there exists no finite CDW order 
that is gapless and incompressible. But for large V the 
CDW order remains finite for incommensurate densities. 
As a result, the region in the phase diagram between 
CDW-I and CDW-II exhibits the coexistence of both the 
diagonal long range order(DLRO) and the off-diagonal 
long range order(ODLRO) which is the signature of the 
supersolid. In order to obtain the boundary that sepa- 
rates the supersolid phase in the phase diagram, we plot 
the Ocdw with respect to V for different densities as 
seen in Fig. [JJ We note that the Ocdw increases sharply 
at some critical value of V, highlighting the transition to 
the CDW phase. To obtain this critical value of V, we 
take the derivatives of Ocdw with respect to V for dif- 
ferent densities in the range 1/2 < p < 1 and p > 1. The 
point where the derivative is a maximum is taken to be 
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FIG. 6: (Color online) Variation of the local density as a func- 
tion of the lattice sites i for doping around p = 1. Panel 
(b) corresponds to the CDW-II phase at p = 1 and panels (a) 
and (c) to the state obtained by removing and adding a boson 
to the CDW-II state. Notice that the crystalline structure is 
preserved for these density changes. V — 9 

the critical point of transition to the CDW phase. The 
order parameters Ocdw as wen as their derivatives as a 
function of V are shown in Fig. [5] The derivative shows 
a negligible peak for p < 1/2, but it shows a sharp peak 
for p > 1/2 indicating the existence of the CDW phase. 
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FIG. 7: (Color online)OcDW as a function p for different V . 
The two peaks at commensurate densities shows the existence 
of CDW-I and the CDW-II phases. The finite order parameter 
for large values of V in the incommensurate density range 
shows the signature of the supersolid phase. 

The phase diagram obtained by plotting the chemi- 
cal potential, p, corresponding to different densities as a 
function of V are given in Fig. 1101 To identify the region 
where the gapped phases exist, we calculate the chemical 
potentials p + and p~ [24], [28| at p = 1/2 and 1 for all 
values of V in the thermodynamic limit and plot them in 
the p — V plane. The boundary of the supersolid phase 
is obtained by calculating the chemical potential p for 
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FIG. 8: The CDW order parameter Ocdw as a function of 
V for different densities in the range 1/2 < p < 1. Note that 
Ocdw increases as V increases. 




V 



FIG. 9: (Color online)CDW order parameters Ocdw (solid 
lines) and their derivatives with respect to V (broken lines) 
are given. For V > Vc ~ 6.4, the derivatives of the order 
parameters show peak at the transition to the CDW phase. 



0.5 < p < 1 and p > 1 at the critical value of V where 
the system enters into the the CDW phase. The critical 
value of V for the transition to the SS phase for density 
close to p — 0.5 is Vc ~ 6.4. This critical value depends 
on density, exhibiting an increase as the system is further 
doped and has a minimum value of Vc = 5.4. For large 
values of V , it is clearly seen that the system continues 
to be in the supersolid phase above CDW-II, while for 
smaller values V , it is in the SF phase. 

The results remain qualitatively similar when U = 5 
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FIG. 10: (Color online)Phase diagram showing all the phases 
for U = 10 
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FIG. 11: (Color online)Phase diagram showing all the phases 
for U = 5 

and the corresponding phase diagram is given in Fig. 1111 
In this case the gapped regions such as CDW-I, CDW-II 
and MI shrink. There is no direct transition from MI 
to CDW-II in sharp contrast to U = 10. Rather there 
are continuous MI-SF and SF-CDW-II transitions. The 
supersolid phase occurs in a small region close to p < 1, 
but the trend is similar to that of U = 10 for p > 1. 



IV. CONCLUSIONS 

In summary, we have obtained the complete phase dia- 
gram for a single species bosonic atoms in the framework 
of the extended Bose-Hubbard model for two different 
values of the onsite interaction U. Our studies have been 
carried out using the FS-DMRG method for a large range 
of densities; 0.25 < p < 1.25. In the large U limit, we 
obtain CDW-I, MI, CDW-II, SF, soliton and the super- 
solid phases and the transitions between them occurring 
at various critical values of the nearest neighbor interac- 
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tion. The supcrsolid phase appears in the density range 
0.5 < p < 1 and p > 1 only in the large V regime. The 
solitons are found to exist for doping above half filling in 
the small V regime and for doping below half filling for 
the entire range of V. For an onsite interaction of inter- 
mediate strength([7 = 5), we find an interesting change 
in the phase diagram. The supersolid phase becomes very 
small in the density range 0.5 < p < 1 and it exists only 
at densities close to 1. 

From an experimental point of view, in addition to the 
optical lattice, there is always a magnetic trap present, 
thereby making these systems inhomogeneous. There- 
fore, it makes it imperative to study this model in the 
presence of a harmonic trap, where all the phases coex- 



ist. Hence it becomes important to look for experimental 
signatures of these phases in the presence of a trap. We 
are currently working in that direction. 
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